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ABSTRACT 

We investigate recent claims that the Hubble sequence of spiral galaxies is scale-free. Funda- 
mental to this investigation is the fact that within the photometric data of 86 spiral galaxies from 
de Jong & van der Kruit (1994) - from which these claims were made - a trend exists between 
morphological type and bulge profile shape. While late-type spiral bulges are described by an 
exponential luminosity profile, the early-type spiral bulges are better described by an r 1 / 2 or 
r 1 / 4 law. Taking the scale-lengths from the best-fitting surface brightness profile models (i.e. 
either using exponential, r 1 / 2 , or r 1 / 4 law profile parameters), we show that in all six passbands 
used (BVRIHK) the early-type spirals have a larger r e /h ratio than late-type spiral galaxies. In 
contrast to this, fitting an exponential profile to the bulges of all spirals results in the mean r e /h 
ratio for the early-type spiral galaxies actually being smaller than the mean r e /h ratio for the 
late- type spiral galaxies (at the 3cr significance level using K-band data). 



Subject headings: 
mation 

1. Introduction 



galaxies: spiral, galaxies: structure, galaxies: fundamental parameters, galaxies: for- 



De Jong (1996b) and Courteau, de Jong, & 
Broeils (1996) have suggested that "the Hubble 
sequence of spirals is scale-free" . They claim that 
"the constant ratio of bulge-to-disk scale-lengths 
appears to be independent of galaxy type" . If 
true, this would not only be at odds with the 
classification scheme posed by Hubble (1926,1936) 
and later Sandage (1961) - in which the bulge-to- 
disk ratio progressively decreases as one goes from 
early to late type spirals (Simien & de Vaucouleurs 
1986) - but would have consequences for theories 
of galaxy formation. 

While the surface brightness profiles of the disks 
of spirals are well described by exponential mod- 
els, the light profiles of the bulges are known to 
possess a range of structural shapes (Andredakis, 
Peletier & Balcells 1995; Carollo, Stiavelli & Mack 
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1998) . These can be easily modelled with the Ser- 
sic (1968) r 1 /™ law. A generalization of de Vau- 
couleurs' (1948) r 1 / 4 law, the free parameter n can 
describe the observed range of bulge profile shapes. 
Indeed, a subset of this model (namely n=l, 2, 
4) was applied by de Jong (1996b) to the surface 
brightness profiles of the bulges in his sample of 86 
face-on spiral galaxies, de Jong found that 60% of 
his sample were better modelled (based on the x 2 
statistic) with an n=l model, while 40% preferred 
a larger value of n, with 15% preferring n=A. Simi- 
larly, Courteau (1996) found, when fitting both an 
n=l and an n—A profile model, that 15% of the 
bulges in his sample of spiral galaxies were bet- 
ter modelled with the n=4 profile (Courteau et al. 
1996). 

Given these results, Courteau et al. (1996) re- 
ported that late-type spirals are best fitted by two 
exponential models, and they chose to represent 
all their spiral galaxies this way. Subsequently, 
their claim for a scale-free Hubble sequence for 
spirals was based upon structural parameters ob- 
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tained from fitting exponential light profile models 
to both the disk and the bulge. 

However, the above percentages become most 
interesting when one notes that the galaxies pre- 
ferring the larger values of n are the early-type 
spirals, while the late-type spirals prefer values of 
n<~l (Andredakis et al. 1995; Moriondo, Giova- 
nardi, & Hunt 1998). Additionally, Andredakis et 
al. (1995) have shown that the bulgc-to-disk ratio 
of luminosities varies systematically with profile 
shape, such that galaxies with larger bulge-to-disk 
luminosity ratio have larger shape parameters. 
Logically, any conclusions drawn from structural 
parameters which have ignored these structural 
differences must surely be questioned (Moriondo 
et al. 1998). By using the best-fitting profile mod- 
els (either n=l, 2, or 4) this paper reinvestigates 
the claim for a scale- free Hubble sequence of spiral 
galaxies. 

2. Data 

We have re-analyzed the data presented by 
Courteau et al. (1996). They presented two data 
sets, however, only one is appropriate for explo- 
rations of galaxy properties as a function of mor- 
phological type. 

Lahav et al. (1995) showed that the dispersion 
in galaxy type index, T, between six experienced 
galaxy classifiers was on average 1.8 T-units, and 
2.2 T-units when comparing the RC3 (de Vau- 
couleurs et al. 1991) T-index with those of the six 
classifiers. A similar figure of disagreement (2.0- 
2.5 T-units) was obtained by four human classi- 
fiers of HST images (Odewahn et al. 1996). Un- 
fortunately, because of this, the larger of the two 
data sets presented in Courteau et al. (1996) - 
243 Sb-Sc galaxies from the 349 Sb-Sc galaxies of 
Courteau (1996) - can not in themselves be used 
to explore possible trends within the Hubble se- 
quence of spiral galaxies. 

What the R-band data of Courteau (1996) does 
enable, is to show that the individual ratios of 
bulgc-to-disk scale-lengths span a broad range of 
values (Courteau et al. 1996, Figure 1). Scale- 
length ratios within just one standard deviation 
of the median are shown to span a range greater 
than a factor of 4, with the lcr confidence interval 
ranging from 0.029 to 0.135, and a long tail in the 
distribution stretching to 0.35. To obtain the ratio 



of the bulge effective radius r e to the disk scale- 
height h, these numbers should be multiplied by 
1.679, giving ratios of r e /h up to ^0.6. Therefore, 
in passing, we stress that caution should be em- 
ployed when using any sort of mean bulge-to-disk 
scale-length ratio, since a broad range of values 
spanning one order of magnitude exists amongst 
the real galaxy population. 

The second data set, that of de Jong & van 
der Kruit (1994), is however useful. It includes 
galaxy types from Sa through to Sm. This sam- 
ple of 86 galaxies actually includes two SO galaxies 
which are removed here as de Jong (1996b) notes 
that their surface brightnesses are significantly be- 
low the trend seen for the rest of the spiral galax- 
ies, and their connection with the early-type spiral 
galaxies is still unclear. The sole irregular galaxy 
(T=10) is also removed, leaving 83 face-on (minor 
over major axis ratios greater than 0.625) Sa to 
Sm galaxies, imaged in six passbands (BVRIHK). 

3. Analysis 

In recent years, some of the limitations of the 
classical surface brightness profile models, such as 
the exponential or the r 1 / 4 law, have been realised. 
Departures in the radial falloff of light from these 
models has been not only detected but success- 
fully modelled for: the dwarf galaxy population 
(Davies et al. 1988; Young & Currie 1994; Binggeli 
& Jerjen 1998), the ellipticals (Caon et al. 1993; 
Graham & Colless 1997), brightest cluster galaxies 
(Graham et al. 1996), and for the bulges of spirals 
(Andredakis et al. 1995). 

The Sersic (1968) law has proved successful in 
parameterizing such departures from the tradi- 
tional models and can be written as 



I{r) = J exp 
= I e exp 



l/n 



(2n- 0.327) 




(1) 



The first line shows how the intensity I varies 
with radius r; 7 is the central intensity where 
r=0. We use hb here to denote the bulge scale- 
length and avoid confusion with the disk scale- 
length which is denoted by h elsewhere in this pa- 
per. The third model parameter, n, describes the 
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level of curvature in the light profile. For exam- 
ple, when n=l the Sersic law is equivalent to an 
exponential light distribution; when n=4 it mim- 
ics the de Vaucouleurs r 1 / 4 law. The value of n is 
of course not restricted to integer values and re- 
mains meaningful up until values of around 10-15. 
The second line is a variant of the first expression, 
with the model parameters now I ei the intensity at 
the radius r e which encloses half of the total light 
of the bulge. Equating like-terms, one has that 
J =/ e exp(2n- 0.327) and {r e /h b )={2n - 0.327)". 
Therefore, when n=l, r e =1.67hb, and when n=2, 
r e = l3.5hi>. One can also easily see why effective 
radii rather than scale-lengths are used for the r 1 / 4 
law, since ft,fc=r e /3466. Given that this paper uses 
parameters from n=l, 2, and 4 Sersic models, we 
have used effective radii rather than scale- lengths. 

de Jong (1996a) fitted three models to the sur- 
face brightness profiles of the bulges, all with an 
accompanying exponential profile model to the 
disk. The goodness-of-fit for each model was mea- 
sured using the \ 2 statistic. 2 For the B, V, H 
and K passbands, it is observed that for every two 
galaxy bulges that are best fit with an n=2 or n=4 
profile, there are three galaxy bulges whose best- 
fitting profile model is the n=l model. For the 
R and I passbands, the number of galaxy bulges 
best fit with the n—1 model equal the number of 
bulges better fitted with the alternative n=2 or 
n=4 models (Table 1). 

In using the best-fitting bulge models (either 
n=l, 2, or 4) the associated model parameters 
were not always reliable. In particular, the r 1 / 4 
model sometimes resulted in values for r e that 
were either inaccurately determined and/or were 
unrealistically large. To accommodate for this, 
each value of r e was inspected and the galaxy 
cither retained, or rejected if Ar e /r e >40% or 
TV Amaa;> 0.5, where r max is the maximum radius 
for which the the surface brightness profiles ex- 
tend (~26±1 in B). This typically resulted in the 
exclusion of only 1 or 2 galaxies from each of the 
morphological type bins T=l-3 and T=7-9 used 
in this comparative study. 

Table 2 shows the difference in the mean value 
of r e /h for the early- and late- type morphological 
class bins used by Courteau et al. (1996). It shows 

2 The data can be found at http://cdsweb.u- 
strasbg.fr/htbin/Cat?J/A+ AS/118/557. 



this ratio for the K- and R-band data fit with 
an exponential bulge model by de Jong (1996b) 
and Courteau et al. (1996). Using the best-fitting 
7i=l, 2, and 4 models, we present this difference 
of means for all six passbands (BVRIHK). How- 
ever, this difference in the ratio is meaningless on 
its own. What is important is the significance of 
this difference, and this depends on the sample 
size and standard deviation of the distributions. 
To this end, we have applied Student's t-test. The 
probability, Prob(t), that the difference in means 
could be as large as it is by chance is given in 
Table 2; small values indicate that the means are 
significantly different from each other. 

4. Discussion 

The majority of the early- type spirals (<Sb) 
prefer to have values of n>l, while late- type galax- 
ies (>Sd) are better fit with an exponential bulge 
(see Table 1). The universal application of the 
exponential fitting function ignores from the start 
real differences in galaxy structure, and introduces 
a systematic bias into the parameterization of 
these galaxies - under-estimating the effective half- 
light radius of the bulge. 3 Figure 1 shows the ratio 
of the effective radii derived from the r 1 / 4 model 
(r e .4) and the effective radii derived from the ex- 
ponential model (r e , exp ), plotted against the ratio 
of the exponential model disk scale length co-fitted 
with the r 1 / 4 bulge model (/14) and the exponen- 
tial disk scale length co-fitted with the exponential 
bulge model (h exp ). It shows that r e ^/r e ^ xp >l, 
while the exponential disk scale-length remains 
largely unchanged as the bulge profile model is 
adjusted. 

Similarly, fitting an n=l profile will over- 
estimate the half-light radii for some of the late- 
type spirals. Although de Jong (1996a) shows 
for the late- type spirals that an n=l model pro- 
vides a better representation of the bulge than 
an n=2 or n=4 model, he also notes that values 
as low as n=0.5 are obtained when applying the 
Sersic profile to the bulge (de Jong 1996a). Fur- 
thermore, Andredakis et al. (1995), in fitting the 
Sersic model to the K-band bulge light profiles of 
30 spiral galaxies, found some Sb-Sd galaxies to 

3 A similar behaviour is known to occur with r 1 / 4 modelling 
of the light profiles of elliptical galaxies (Graham & Colless 
1997, Figure 11). 
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have bulge profiles with shape parameters smaller 
than 1. Consequently, restricting the structural 
profiles of all late-type galaxies to be described by 
an n—1 model may be increasing their mean bulge 
scale-length and hence reducing the true difference 
between the r e /h ratio of the early- and late- type 
spirals. That is, the probabilities in Table 2 may 
be larger than they should. 

As stated by de Jong (1996b), K-band data is 
the passband of choice for such studies, making 
it "possible for the first time to trace fundamen- 
tal parameters related to the luminous mass while 
hardly being hampered by the effects of dust and 
stellar populations." Indeed, some of the galaxies 
in de Jong sample were noted to possess dust lanes 
and circumnuclear star formation. Furthermore, 
bulges are brighter in K than in B with respect 
to the disk, and since the bulge/disk decomposi- 
tion is easier when the bulge is relatively brighter, 
the fitting algorithm therefore works better in the 
K-band (de Jong 1996b). 

Using exponential bulge models, Courteau et al. 
(1996) mention that the r e /h ratios of the early- 
type spirals appear systematically below the aver- 
age r e /h value for all spiral galaxy types. They as- 
sert that this difference is not large, and claim that 
the constant ratio of bulgc-to-disk scale-length is 
independent of Hubble type. However, our analy- 
sis of the exponential models fitted to the K band 
data of de Jong (1996b, Figure 18) reveals that the 
mean value of r e /h for the Sa-Sb type galaxies is 
actually smaller than that for the late-type spirals 
at a significance of 98% (3cr)! (Table 2). Similarly, 
with the R-band data presented by Courteau et al. 
(1996), and in fact for all wavelengths used (ex- 
cluding the V-band), the ratio of r e /h is smaller 
for the Sa-Sb galaxies than it is for types >Sbc. 
This result is at odd with the classical picture 
of the Hubble sequence, where early-type spirals 
have larger bulge-to-disk scale-length ratios than 
late-type spirals. 

Due to the use of exponential bulge models 
for the Sa-Sb type galaxies, the above result can 
be understood in terms of systematically under- 
estimating the size of these bulges. Correcting for 
this, by taking the best-fitting structural param- 
eters, from either the n—1, 2, or 4 models, we 
find that the situation reverses itself. The aver- 
age value of r e /h for the Sa-Sb type galaxies is 
found to be larger than the average value of r e /h 



for galaxy types >Sd, in all six passbands. Ta- 
ble 2 shows that the probability that the Sa-Sb 
type galaxies have the same mean r e /has the Sd- 
Sm type galaxies is weakly ruled out at the 1.5-2 
o level in five of the six passbands used by de 
Jong & van der Kruit (1994). Interestingly, it is 
the K-band data which suggest that the difference 
in means is not significant. However, this result 
in itself is significant when compared to the re- 
sult obtained using only exponential bulge profile 
models. Using the best-fitting models, the aver- 
age r e /h ratio for the sample is larger - at the 3 o 
significance level - than when obtained using only 
the n=l model. 

We plan to refine this work by fitting a Ser- 
sic profile with free (i.e. not fixed) shape param- 
eter, n, to the bulges of the spirals in the sample 
of de Jong & van der Kruit (1994). Furthermore, 
Courteau et al. (1996) noted that about 1/3 of the 
sample of galaxies from de Jong (1996a) had a bar 
modelled as an additional component - requiring 
eight structural model parameters for these galax- 
ies. While de Jong modelled a bar when fitting the 
exponential bulge models to the 2D images, his 
one-dimensional decomposition technique which 
he used to fit the r 1 / 4 and r 1 / 2 bulge models did 
not allow for the influence of a bar. Subsequently, 
we must caution that failure to model the bar in 
the ID data used here may influence the scale- 
lengths obtained. 

Arguments for secular evolution, namely the 
exponential bulge light profile and the restricted 
range of bulge-to-disk scale-lengths, are either 
wrong or questionable. Andredakis et al.'s (1995) 
alternative to secular evolution - based upon the 
continuous trend between galaxy structure, as 
measured by n, and galaxy type - is largely sup- 
ported by the data of Courteau et al. (1996). In 
the framework of this model, n-body simulations 
(Andredakis 1998) have shown how the imprint of 
disk formation is left upon the bulge, creating the 
observed trend between shape parameter and mor- 
phological type. Yet another alternative is offered 
by Aguerri & Balcells (1999), where the shape of 
the bulge grows from an n—1 profile to larger val- 
ues of n as shown through n-body simulations of 
merger events. 

Whether the bulges of spiral galaxies formed 
after the disk, as in the secular evolution model 
(Courteau et al. 1996), or, whether the bulge is 
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in fact older than the disk (Andredakis 1998, and 
references within) may be better answered when 
the range and trends of bulge-to-disk ratios are 
better known. 

We thank Marc Balcells for his comments on 
this paper prior to its submission. We also wish to 
thank the anonymous referee for their comments 
and suggestions. 
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Fig. 1 . — For each galaxy, the value of r e from the 
n=4 bulge model, divided by the value of r e from 
the n=l bulge model, is plotted against the corre- 
sponding ratio of h from the two associated expo- 
nential disk models. While the scale-length of the 
exponential disk remains largely unchanged, the 
effective radius of the bulge is shown to be larger 
with an r 1 / 4 bulge model than with an exponential 
bulge model. (K-band data.) 



Table 1 
Best-fitting bulge models 



Morphological 


No. of galaxies best fit with each model 


Type 


Exponential ( 


r l/2 )r l/4) 


(no data) 




B-band 






Od,iJd,U,ijlJ 


to 


12 





uDC,ijC,ijC(1 


31 


16 


1 


SHtti Sim 


11 


2 







V-band 








7 


8 


7 


u Ut,iJt,OtQ 


24 


18 


6 


Sd,Sdm,Sm 


10 


2 


1 




R-band 






Sa,Sab,Sb 


6 


16 





Sbc,Sc,Scd 


25 


21 


2 


Sd,Sdm,Sm 


10 


3 







I-band 






Sa,Sab,Sb 


5 


11 


6 


Sbc,Sc,Scd 


21 


24 


3 


Sd,Sdm,Sm 


10 


3 







H-band 






Sa,Sab,Sb 


6 


8 


8 


Sbc,Sc,Scd 


13 


5 


30 


Sd,Sdm,Sm 


4 


2 


7 




K-band 






Sa,Sab,Sb 


8 


13 


1 


Sbc,Sc,Scd 


27 


18 


3 


Sd,Sdm,Sm 


12 


1 






Note. — This table lists the number of spiral galaxy 
bulges whose surface brightness profiles are best fit with 
an exponential model, and the number of galaxy bulges 
that are better modelled with the alternative (r 1//2 or 
r 1 / 4 ) profile. The fourth column displays the number 
of galaxies for which no \ 2 goodness-of-fit estimate was 
available, \ 2 data taken from de Jong (http://cdsweb.u- 
strasbg.fr/htbin/Cat?J/A+AS/118/557). 
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Table 2 

Comparison of the r e /h data distributions 



B ^d <lf) 1 -<lf) 2 Prob(t) 

Exponential bulge model data 
l:(SO,Sa,Sab,Sb) vs. 2:(Sbc,Sc,...Sm,Irr) 
K-band (0.119-0. 162)=-0.043 02% 
R-band (0.112-0.124)=-0.012 34% 

Exponential bulge model data 
l:(Sa,Sab,Sb) vs. 2:(Sd,Sdm,Sm) 
K-band (0.118-0.195)=-0.077 03% 
R-band (0.111-0.126)==-0.015 52% 

Best-fitting bulge models 



l:(Sa,Sab,Sb) vs. 2:(Sd,Sdm,Sm) 



K-band 


(0.240-0.228)= 


=+0.012 


79% 


H-band 


(0.485-0.173)= 


=+0.312 


07% 


I-band 


(0.348-0.190)= 


=+0.158 


24% 


R-band 


(0.333-0.186)= 


=+0.147 


23% 


V-band 


(0.671-0.161)= 


=+0.510 


14% 


B-band 


(0.321-0.160)= 


=+0.161 


18% 



Note. — Comparison of the r e /h data dis- 
tributions for different morphological-type 
bins. Column 1 shows the passband used. The 
difference between the mean values from the 
two distributions (as listed 1: and 2: in the 
table sub-headings) is shown in column 2. Col- 
umn 3 gives the significance that the two dis- 
tributions have the same mean value, as de- 
rived from Student's t-test and allowing for 
different population variances between the two 
data sets. Small probabilities indicate that the 
data sets are different. 



8 



T 1 1 1 r 



a 



CD 



6 



_i i i i_ 



_i i i i_ 



ID 



